Two new topologically ordered glass phases of smectics confined in anisotropic 

random media 
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Smectics in strained aerogel exhibit two new glassy phases. The strain both ensures the stability 
of these phases and determines their nature. One type of strain induces an "XY Bragg glass" , while 
the other creates a novel, triaxially anisotropic "m=l Bragg glass". The latter exhibits anomalous 
elasticity, characterized by exponents that we calculate to high precision. We predict the phase 
diagram for the system, and numerous other experimentally observable scaling laws. 
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Liquid crystals confined in random porous structures, 
have become a subject of considerable interest. (!]] A re- 
cent theoretical study unambiguously demonstrated that 
conventional (quasi-) long-ranged smectic order is im- 
possible in 3d in the presence of (even arbitrarily weak) 
quenched pinning imposed by these random structures, 
e.g., aerogel. Q It was proposed that a positionally dis- 
ordered but topologically ordered "smectic Bragg glass" 
(SBG) phase would become the new thermodynamically 
distinct low-temperature phase in these smectic systems. 
However, for quenched random isotropic structures it was 
impossible to make a compelling theoretical argument for 
the stability of such a glass phase. 

In this Letter, we show that we can make such 
a compelling argument for smectics in a uniaxially 
strained aerogel, which certainly exhibit two types of 
low-T BG phases, that are thermodynamically distinct 
from the high-T nematic (or perhaps "nematic elastic 
glass" (NEG)) and isotropic liquid phases. For paral- 
lel nematogen-surface alignment (assumed throughout), 
a stretch (Fig.[l]a) of the aerogel will lead to an "XY-BG" 
in the isotropic universality class of randomly pinned 
vortex lattices, CDW's, and random field XY magnets 
(RFXY), [H while a compression (Fig.[l]b) will lead to a 
novel "m = 1 BG" , with triaxially anisotropic scaling, 
that should be similar to that of a discotic in isotropic 
aerogel [Q. For homeotropic alignment, the phases re- 
verse with respect to stretch and compression. 
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FIG. 1. (a) Stretch along the z direction, (b) Compression 
along _L direction. 

We predict two possible low, constant-T phase dia- 
grams, depending on whether the SBG is stable (Fig.Ep) 



or not (Figj^a). Recent experiments || suggest the for- 
mer possibility. The loeii of the phase boundaries in 
Fig.Ha, for small strain, a, are universal and satisfy 



A(cr) cx (K 3 B)^ 2 {a/B) p , 



(1) 



where a is proportional to the uniaxial stress applied to 
the aerogel fibers, A is a measure of the tilt disorder, B 
and K are bulk smectic elastic moduli, and p a univer- 
sal exponent expressible in terms of anomalous elasticity 
exponents f\B and t)k for unstrained aerogel. Our best 
estimate is p m 1/3 in 3d. || 

Our model of the smectic in aerogel treats the local 
smectic layer displacement u(r) and the local nematic di- 
rector n(r) as the only important fluctuating quantities, 
ignoring fluctuations in the magnitude of the smectic 
order parameter ip — \ip\e lq ° u ^ about its mean \ip \- 
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FIG. 2. Two possible phase diagram topologies, depending 
on whether S-BG is stable for isotropic confinement (a = 0). 

The important effects of the aerogel are com- 
pletely described by only two disorder types. One 
is the random-field translational disorder SH r f = 
Tie $d d r\ip \V{v)e lq ° u{ - r \ coupling to u(r), with V(r) a 
complex random potential which at long scales can be 
accurately represented as zero-mean and short-ranged 
Gaussian statistics with V{r)V*(r') = Ay<5 d (r-r'). The 
other type of disorder is the random-tilt orientational 
disorder given by SH t = — fd d r(g(r) ■ n) 2 , describing 
the tendency of nematogens n(r) to align along the lo- 
cal aerogel strand directed along g(r), and at long-scales 
is completely described by short-ranged [^) correlations 
gi(r)gj(r') — l/2(^/~K5ij— jeiej)S d (r— r'), with e the uni- 
axial direction (i.e., the axis of the strain applied to the 
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aerogel). In the above Ay = T u (a f /L f ) d ~ dF (\/(a f q )), 
A = T n (a,f /Lf) d ~ dF , df is the aerogel's fractal dimen- 
sion for scales a/ < r < Lf, and T u , T n , 7 are phe- 
nomenological parameters, with 7 the anisotropy param- 
eter which at small strains is proportional to yA and 
the stress a applied to the aerogel. 7 < for a stretch 
illustrated in Fig.[l]a. 

Assuming (as we'll verify a posteriori) that fluctuations 
in n from a perfect alignment with the smectic layer nor- 
mal (taken along z) are small, allows us to integrate n out 
of the partition function, with the only effect of replacing 
8n — > Vj_u. |2| The resulting Hamiltonian is given by 



H 



f(Vi U ) 2 + f(^-i(V^) 2 ) 2 
+ (<? z (r)Vyu) 2 - (g(r) • V ± u) 2 - 2 52 (r)g(r) • V ± u 



\^ \TZe{v(r)e^ u ^} 



(2) 



The form of the anharmonic elastic terms is dictated by 
the underlying invariance of the bulk smectic phase under 
rotations about any axis lying in the ry -plane. 

After introducing n replica fields u a and integrating 
out the disorder, we obtain the Hamiltonian whose form 
strongly depends on the type of uniaxial strain. Stretch- 
ing the aerogel strands will cause the layer normal, z, 
to align with e (Fig.[j]a). Smectics confined inside this 
structure, to harmonic order in elasticity (with elastic 
anharmonicity irrelevant) are described by 



H 



XY 



L Jr a =l 
9T / 



K(Vlu a ) 2 + B(d z u a ) 2 + \-f\(V ± u a ) 2 

(A + VA| 7 |)(Vx«a)-(V_Lti i8 ) 



r a,f3=l 
+ Ay cos[q a (u a - Up)] 



(3) 



where Ay = 2|^ | 2 Ay. At scales smaller than a 
crossover scale £j_ (to be defined below), the behavior is 
that of a smectic pinned by isotropic unstrained aerogel. 
|^[ On longer scales, however, such anisotropically pinned 
smectic crosses over to isotropic scaling behavior of the 
RFXY model. We therefore predict that smectics pinned 
by such anisotropic weak disorder will exhibit the XY-BG 
phase, with its universal disorder-induced logarithmic 
layer wandering character, (u 2 (r)) = C (d) (in r) / q 2 . jj| 
However, unlike 3d bulk smectics, which show the famous 
Landau-Peierls thermally-dxivcn lnr fluctuations, here 
C(d) is universal, the logarithm persists in all 2 < d < 4, 
and smectic layers are pinned. The immediate conse- 
quence is that X-ray scattering will exhibit real-space 
power-law decay (pg(t)p-c(0)) oc r _ '^ G ' with a univer- 
sal 77(G) exponent (G = mq ). 

If, instead, the aerogel is uniaxially compressed, i.e., 
7 > 0, we expect that one of the (previously soft) rj_ 



smectic axes (x or y) will orient along the axis of com- 
pression eh (Fig.|l|b). We denote this e/j-directed axis as 
"hard' (h), and call the other _L axis, orthogonal to e^, 
"sq/T (s) axis, i.e., ry = (r^,r s ). The resulting effective 
Hamiltonian describing this system at long scales is 



H m=1 = ~J K (viu a ) 2 + B(d z u a - hy^u a ) 



A 



«{V h u a ) 2 - y(^±w Q )-(V 



a, 0=1 



up) 



(4) 



where we have neglected the positional random-field dis- 
order, Ay, which can be shown to be subdominant at 
long length scales. Q H m= i, Eq.[|, implies that the non- 
interacting propagator G aj g(q) = V -1 (u a (q)up(~ q))o = 
TG(q)S aP + Aq 2 ± G(q) 2 , with G(q) = l/{Kq\ + iq \ + 
Bq 2 ). As usual, at long length scales, the disorder (A) 
contribution to layer roughness dominates over the ther- 
mal (T) part of G a p (q) . 

We first note that for vanishing strain 7 > 0, or equiv- 
alently at very short length scales, H m —\ and the cor- 
responding propagator reduce to those characterizing a 
smectic in unstrained isotropic aerogel. Q The asymp- 
totic long scale behavior of the full model described by 
H m= i is reached via two independent crossovers from 
the Gaussian, unstrained fixed point, during which the 
aerogel anisotropy 7, and the nonlinear elasticity, re- 
spectively become important. The qualitative form of 
this crossover is determined by the relative magnitudes 
of the corresponding bare couplings. For sufficiently 
weak strain (7 < 7c ), the elastic anharmonicity be- 
comes important first and this occurs at a crossover 

length scale £± L oc ( b'i/^aX ^ 1 determined by the smec- 
tic in unstrained aerogel. [0 In this case the system first 
crosses over from the unstrained Gaussian to unstrained 
anomalous fixed point. The final crossover to asymp- 
totic strained anomalous behavior takes place within the 
anomalously elastic smectic described by the wavevector- 
dependent elastic constants j| and occurs at qy such 
that -K"(qy)gy ~ 7<7y, with if (qy) calculated in Ref. [Q, 
i.e. at CI ~ [A7(7(£^ i )) f >*-] 1 /(2->)ic) ( we use tilde sym- 
bol for exponents for isotropic disorder). 

For the remainder of this paper we will focus on the 
other crossover scenario in which the strain 7 is suffi- 
ciently large (7 > 7c ) that the crossover from Gaussian 
unstrained to Gaussian strained elasticity takes place at 
£y = y/ A/7, before elastic nonlinearities become impor- 
tant. The critical value of 7 that delineates between these 
two crossover scenarios is 7 C = K/{£± L ) 2 . 

For 7 > 7c, on scales longer than = yKpy the ef- 
fective Hamiltonian (and the propagator G derived from 
it) is identical to that given in Eq.|], but with all Vy re- 
placed by V s , with r s a subset of ry axes remaining soft 
even in the presence of aerogel anisotropy. Our goal then 
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is to assess the role of elastic nonlinearities, at this new 
strained Gaussian fixed point, which become important 
beyond an even longer nonlinear crossover length scale 
(along the 'soft' direction). 
The scale £^ can be determined from a simple per- 
turbation theory in these nonlinear couplings of H m —i, 
and is the length at which the effects of anharmonic elas- 
tic terms become significant. For example, the diagram- 
matic correction to the bulk modulus B, due to these 
elastic nonlinearities, is given by 



SB(L) 



2 

2tt(7 - 2d) 



[TG(q) 2 



2Ag 2 G(q) 3 ] qt , 



B :i 



d-2j(7-d 



1/2 



L 



7-2d 



(5a) 



(5b) 



where in the above we have kept only the dominant dis- 
order infrared divergent contribution, cutoff these long 
scale divergences by q s > 1/L, and analytically con- 
tinued to arbitrary dimension d, with a single smec- 
tic ordering coordinate z, a single soft coordinate r s , 
and d — 2 hard axes with coordinate r/,. The constant 
C d = 2ir d / 2 / ((2ir) d T(d/2)) and d = r(d/2)r(3-d/2)/2. 
For d < d uc — 7/2, the corrections to B, Eq.[3l], grow 
with cutoff L and become significant for scales L > 
such that \SB(^ L )\ = B, signaling the breakdown of 
conventional harmonic elasticity. We find 



NL 



2^(7-2rf)if( 7 - d >/ 2 7 ( d - 2 )/ 2 ^ 1/(7 - 2d) 



(6) 



The corresponding lengths along the z and h axes are 
given by g L = (Z? L ) 2 /\ B and ^ L = (£f L ) 2 /A 7 , where 
X B = {K/B) 1 / 2 and A 7 = (if/7) 1/2 . Identical crossover 
lengths scales are obtained if one instead studies pertur- 
bative corrections to K or A. 

To go beyond these crossover length scales we use 
the renormalization group (RG), which consists of inte- 
grating out short-scale modes, perturbatively in elastic 
nonlinearities, and rescaling the lengths and long wave- 
length part of the fields with r s = r' s e e , = r' h e u ' h , 
z = z'e^" and u a (r) — e xl u' a (r'), so as to restore the 
uv cutoff back to A ~ 1/a. The underlying rotational in- 
variance ensures that the graphical corrections preserve 
the rotationally invariant operator (d z u a — ^(V s it a ) 2 ), 
renormalizing it as a whole. It is therefore convenient 
(but not necessary) to choose the dimensional rescaling 
that also preserves this operator; the appropriate choice 
is x = 2 — oj z . 

Using the above defined analytical continuation in d, 
RG to one-loop order, gives the following flow equations 



dB{t) 

d£ 
dK(£) 

dt 



= f<j + (d- 2)uj h - 3lo z 
= (1 + (d- 2)u) h - lo z - 



3ff \ 

9 
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8\/2' 



K 



(7a) 
(7b) 



d(A/T)(£) 
dl 



= (3 + (d - 2)u h -u z + -A=) (A/T) , (7c) 



where we've defined a dimensionless measure of disorder 
g = A(B/(K 7 - d Y~ 2 ) 1/2c d-ik 2d - 7 , which flows accord- 
ing to 



dg(i) 
dl 



2eg 



15 
64^ 
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(8) 



with e = 7/2 — d. Because all relevant anharmonic terms 
in H m= i appear with V s , there are no graphs correcting 
7 and therefore no anomalous 7 elasticity to all orders. 
As required, the flow of g is independent of the arbitrary 
choice of the anisotropy rescaling exponents u>h and w z . 
The growth of g for d < d uc = 7/2 is an indication that 
the long scale properties of our system, even at a finite 
temperature T, are dominated by disorder. The even- 
tual termination of this flow at a nontrivial, glassy T = 
fixed point g* = el28\/2/15, leads to strong disorder- 
generated power-law anomalous elasticity. 

One consequence of the anomalous elasticity is that the 
long-scale elastic constants K, B, and disorder variance 
A become wavevector dependent. 



K(k) = Kk~^f K (k h /k<»,k z /k$- 
B(k)=Bkrf B (k h /k^,k z /k^) 
A(k) = Afc s -"-/ A (k h /k^,k z /k^) 



(9a) 
(9b) 
(9c) 



7(k) = 7 , with the anisotropy exponents C, z = 2 — (r] B ■ 
T]k)/2 and £h = 2 — Vk/2. The exponents obey 



t\b 7 — d 
7 - 2d + r) A = — H Tjjc 



(10) 



exactly, due to the underlying exact rotational invari- 
ance of Eq.^ about e^. To leading order in e = 7/2 — d, 
m = g*/8V2 = 16e/15 = 8/15, r/ B = 3g* /32V2 = 
12e/15 = 2/5, and r)A = g*/32V2 = 2e/15 = 2/15, the 
last equalities holding in d = 3 (e = 1/2). Since e = 1/2 
is quite small, we expect these exponents to be quantita- 
tively accurate. 

The RG e = 7/2 — d expansion treatment presented 
above is nicely complemented by an e = 4 — d expan- 
sion arising from a different analytical continuation to d 
dimensions, in which there are d — 2 soft coordinates r s 
and only a single hard axis. The corresponding expo- 
nents are given by t)k = 3e/8 = 3/8, f\ B — 3e/4 = 3/4, 
and f] a = e/8 = 1/8, with good agreement in d = 3 
(except for rj B ) with the e = 7/2 — d expansion results. 
The exact exponent relation for the e = 4 — d expansion 
is given by 4 — d + ??a = Vb/2 + 2f]K and reassuringly 
agrees with Eq.[l0] in d = 3. 

Further accuracy can be gained by weighted-averaging 
of the 7/2 — d and 4 — d expansions, according to: 
t]k,a — * (4?7_ftr,A + Vk a)/^- The factor of 4 reflects the 
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higher accuracy of the 7/2 — d = e expansion. The predic- 
tion for rjB is then made using the exact scaling relation, 
giving, in d = 3, t)k = 0.50, r/s — 0.26, and ?7a = 0.13. 

We now study the translational and orientational order 
in the presence of this strong disorder-driven anomalous 
elasticity. The former is characterized by the growth of 
smectic layer roughness with e.g., r s 

C(r s ) = {(u(r s ,0 h ,0 z )-u(0 s ,0 h ,0 z )) 2 ) , 

/d 3 k 
— 2(1 - cos(A ;s r s ))A(k)fciG 2 (k) , (11) 

from which the translational correlation length £^ (the 
inverse of the X-ray diffraction peak width) , can be com- 
puted via the condition C(r s = £^ ) = a 2 , with a the 
smectic layer spacing. £y is determined by the relative 
order of many crossover length scales. For < £^ 
is identical to that due to isotropic disorder. Q In the 
opposite regime, > there are three possibilities, 
depending on whether the anomalous elasticity sets in be- 
fore or after the layer roughenness reaches a, and whether 
the isotropic-to-anisotropic crossover takes place near the 
harmonic or the anomalous elastic fixed point: 

-e L (^) 2 , £i<e<ef L ,7>7c 

,x = I ^ L (^)^f^> G < ^ L < 7 > 7c 

s | f c fa\ {VB 1 VK) x 

t (^ , £r < & < , 7 < 7c 

(12) 

From these correlation lengths we see that it is the ratio 
a j Xb which determines whether £f lies in a length scale 
regime in which anharmonic effects are important. For 
small B, Xb S> a and anharmonic effects are unimpor- 
tant. Note also that in the strained length scale regime, 
£s will depend on B,K,A and 7. Thus, one could test 
the predictions of Eqjl^ by measuring the dependence of 
£g on the strength of compression (i.e., 7) which could 
be adjusted directly. In all length scale regimes, the X- 
ray correlation length is finite even as T — > 0, signaling 
the destruction of the conventional (quasi-) long-ranged 
translational smectic order. 

As emphasized in Ref. ||, this lack of translational 
order does not imply that the low temperature phase 
replacing the smectic is simply a nematic (or "worse" , 
isotropic). Our detailed calculations Q indicate that, 
in fact, despite the lack of the (quasi-) long-ranged 
smectic order dislocation loops remain bound for weak 
anisotropic disorder, and therefore the low temperature 
phase replacing the smectic must be distinct from the ne- 
matic, separated from it by a thermodynamically sharp 
dislocation unbinding phase transition. We call this pu- 
tative low-temperature phase the "m = 1 -Bragg Glass" . 

The stability of this exotic glass phase is contingent 
upon our implicit assumption of long-ranged orienta- 
tional (nematic) order. That this assumption is valid can 



be easily seen by computing (|<5n| 2 ) = (|Vu| 2 ), and tak- 
ing into account the wavevector dependent elastic mod- 
uli K (k) and B(k) and disorder variance A(k), as given 
by Eqs.[|. There are unstrained and strained contribu- 
tions to (|<5n| 2 ), arising from modes with q± > 
and q± < respectively. Using the corresponding 

anomalous exponents tjk, Vb, Va, |§ in the computation 
of the strained and unstrained parts, for finite strain 7 
and weak disorder A, we indeed find long-ranged orien- 
tational order. In the weak strain limit (7 < 7 C ), the 
unstrained part dominates, in 3d growing in a universal 
way with decreasing strain 7 and increasing disorder A 
as (A^/j^ 1 ), where n — rj B /{2 — f] K ). Using Ref. |§] we 
estimate /i to be 3/2. Since ( | <5n| 2 ) can therefore get ar- 
bitrarily large at small 7 and large A, we expect our sys- 
tem to be in the orientationally disordered liquid phase in 
this range of parameters. On the other hand for large 7 
and small A the system will exhibit long-ranged orienta- 
tional order and as illustrated in Fig.||a will therefore be 
in the m — 1 BG phase. In analogy with the Lindcmann 
criterion for melting, the phase boundary is roughly de- 
termined by the condition (|<5n| 2 ) w 0(1). This leads to 
the phase boundary quoted in Eq.|l|, and illustrated in 
Fig.|a. 

On the other hand, rather than rely on the untrustwor- 
thy (in 3d) 5 — e expansion, which predicts no SBG for 
isotropic (a — 0) disorder 0, we can infer the topology 
of the phase diagram based on the preliminary experi- 
mental evidence H, which suggests the stability of SBG 
for weak isotropic disorder. This suggests that the m = 1 
BG extends all the way down to vanishing strain, a > 0, 
as illustrated in Fig.^>. 

Light scattering, which measures director correlation 
(Sm(q)Snj(— q)) provides an independent means to test 
the predictions of the theory. Finally, since anomalous 
elasticity also implies a nonlinear stress-strain relation at 
arbitrarily weak stress, our predictions for it can be inde- 
pendently probed in an a.c. acoustic experiment, search- 
ing for an unusually large second harmonic response. 
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